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The computational cost of numerical simulations of QCD with 
light dynamical Wilson-quarks is estimated. The qualitative be- 
^ ! haviour of the pion mass and coupling at small quark masses is 

discussed. 
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: 1 Introduction 

The eightfold way of hadron spectra is the consequence of the SU(3) flavour 
symmetry in QCD which follows from the existence of three light quarks 
{u, d and s). Numerical simulations of QCD have to take into account this 
basic fact. The "quenched approximation", which corresponds to the limit 
of infinitely heavy virtual quarks, is not reliable because of its unknown 
systematic errors (see, for instance, 0]). 

Chiral perturbation theory (ChPT) is very useful for extrapolating 
the results of numerical simulations to the small u- and (i-quark masses but, 
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of course, it has a finite range of applicability. In practice this means that 
one has to stay in a region where the one-loop ChPT formulas give a good 
approximation. For investigating the masses and couplings of pseudoscalar 
mesons one should reach the quark mass range below about one quarter of the 
strange quark mass {mud < 1"?.^) @, Present large scale QCD simulations 
- especially with Wilson-type quark actions - are not yet at small enough 
quark masses for the applicability of chiral perturbation theory. This is 
clearly displayed, for instance, by the absence of curvature ( "no chiral logs" ) 
in the present data about the dependence of and f-,, versus the quark 
mass rriq (see, for instance, 0, 0) B 0)- 

Going to light quark masses in unquenched QCD simulations is a great 
challenge for computations because known algorithms have a substantial 
slowing down towards small quark masses. The present status has been 
summarized by the contributors to the panel discussion at the Berlin lattice 
conference |T^. In these contributions the computational costs were estab- 
lished at relatively large quark masses and extrapolated into the interesting 
range of small quark masses. Recent developments show, however, that at 
small quark masses new difficulties appear which might invalidate this extrap- 
olation |Tl|, |T2|. This was the motivation of our collaboration to investigate 
in a series of algorithmic studies the computational cost of numerical sim- 
ulations down to about one fifth of the strange quark mass {rriud < k^s) 



13, 14|. In our tests we used the two-step multi-boson algorithm |15]. 



2 Cost estimates 

The computational cost of obtaining a new, independent gauge configuration 
in an updating sequence with dynamical quarks can be parametrized, for 
instance, as 

C = F (^)"(^)-. (1) 

Here ro is the Sommer scale parameter, the pion mass, L the lattice 
extension and a the lattice spacing. The powers -Zyr.L.a and the overall constant 
F are empirically determined. 

In order to limit the necessary computer time our collaboration first per- 
formed high statistics simulations with Nf = 2 degenerate quark fiavours on 
coarse lattices (8^ • 16) with a lattice spacing a ~ 0.27 fm. (a is determined 
by To — 0.5 fm from the mesured value of ro/a.) This gives for the space- 
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Quark mass dependence of plaquette autocorrelations 
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c= 7.92(68) 
z = -2.02(10) 
X^/ndf = 1 .8 
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^ M, = (rom/ 



Figure 1: Power fit of the plaquette autocorrelation given in units of 10^ • 
MVM as a function of the dimensionless quark mass parameter Mj.. The 
best fit of the form cM^ is at c = 7.92(68), z = -2.02(10). 



extension L ~ 2.2 fm, which is sufficient for keeping finite volume effects 
tolerable, for instance, for and fj^. For the definition of the quark mass 
the dimensionless quantity 

Mr = (rom,)2 (2) 

has been used, which already appears in (|I]). In this variable the strange 
quark mass can be defined as Mr, strange = 3.1. 

Our results for the integrated autocorrelations of different physical quan- 



tities (Tint) have been described in detail in |T^. For illustration see figures 
|I| and where the quark mass dependence of t-pI^"?"^**^ and r™'^'' are shown, 
together with power fits. The units for the costs are 10^ fermion-matrix- 
vector- multiplications (MVMs). Comparing the two figures one sees that the 
"cost" for is substantially less than the one for the average plaquette. 
The fitted powers of the quark mass are also different. In the notation of the 
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Quark mass dependence of pion mass autocorrelations 



fit:c*x^ 

c= 1.99(16) 
z = -1.47(1 6) 




Figure 2: Power fit of the autocorrelation of the pion mass given in units 
of 10^ • MVM as a function of the dimensionless quark mass parameter M^. 
The best fit of the form cM^ is ate = 1.99(16), z = -1.47(16). 



formula in (|I|) gives z^, ~ 4, whereas r^^" is consistent with Zj, ~ 3. 

The pion couphng constant /^r, which is related to the constant in front of 
the pion contribution in the pseudoscalar propagator, has even shorter au- 
tocorrelations than mTT- This is obviously rather advantageous for obtaining 
physical information from the quark mass dependence of the pion mass and 
coupling (see |^). 

Our collaboration also studied the volume dependence of the cost by 
choosing a few 8^-16 runs and repeating them on either 16^ or 12'^ -24 lattices. 



without changing other parameters ||1J] . The parameters of the chosen points 



are collected in table |T] and some results are shown in table 0. The error 
analysis and integrated autocorrelations in the runs on the larger lattices have 
been obtained using the linearization method of the ALPHA collaboration 
| T6| . The conclusion from these tests is that the cost increase with the lattice 
volume is quite acceptable because it is close to the trivial volume factor. In 
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Table 1: Runs for comparing the simulations costs at different volumes. For 
the definition of algorithmic parameters see 



run 


lattice 


P 


K 


ni 


^2 


^3 


A 


e- 10^ 


(e) 


8^ • 16 


4.76 


0.190 


44 


360 


380 


3.6 


2.7 


(el6) 


16^ 


4.76 


0.190 


72 


350 


440 


3.6 


2.7 


(f) 


8^ ■ 16 


4.80 


0.190 


44 


360 


380 


3.6 


2.7 


(fl2) 


123 . 24 


4.80 


0.190 


72 


500 


560 


3.4 


1.36 


(h) 


8^ • 16 


4.68 


0.195 


66 


900 


1200 


3.6 


0.36 


(hl6) 


16^ 


4.68 


0.195 


96 


860 


1100 


3.6 


0.36 



case of the autocorrelation of the pion mass the observed increase turns out 
to be even smaller. This is partly due to the intrinsic fluctuation present 
in the pion propagator which is originated from the freedom of randomly 
choosing the position of the source. 



3 Chiral logs? 

The behaviour of physical quantities, as for instance the pseudoscalar meson 
("pion") mass m^r or pseudoscalar decay constant /^r as a function of the 
quark mass are characterized by the appearance of chiral logarithms. These 
chiral logs, which are due to virtual pseudoscalar meson loops, have a non- 
analytic behaviour near zero quark mass of a generic form mglogmg. They 
imply relatively fast changes of certain quantities near zero quark mass which 
are not seen in present data P, [7|, p|, 



In our range of quark masses {nig ^ the logarithmic curvature 

due to chiral logs has to appear, at least in the continuum limit. Although 
we have rather coarse lattices (a ~ 0.27 fm) and, in addition, up to now 
we are working with unrenormalized quantities - without the Z-factors of 
multiplicative renormalization - it is interesting to see whether the effects of 
chiral logs can already be seen in our data. 

Let us recall the one-loop ChPT formulas for m^: 
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Table 2: Results from the runs with parameters shown in table [|. The lat- 
tice extension L is obtained from the value of ro/a. is the quark mass 
parameter defined in The cost of an update cycle Cue one? the integrated 
autocorrelations of the average plaquette and of the pion mass are given, 
respectively, in 10^ MVMs and flops. 



run 


L[fm] 


Mr 


Cue [kMVM] 


[flop] 


O [flop] 


(e) 


2.31(6) 


1.473(88) 


8.50 


4.59(37) -10^^ 


1.94(31) -10^^ 


(el6) 


4.57(9) 


1.401(79) 


12.4 


7.5(1.3) -10^^ 


5.02(55) -10^3 


(f) 


2.25(4) 


1.026(51) 


8.48 


7.47(84) -10^3 


1.76(59) -lO^^^ 


(fl2) 


3.02(9) 


1.37(9) 


12.3 


2.40(41) -lO^^ 


4.52(82) -lO^^ 


(h) 


2.27(5) 


0.806(50) 


16.2 


1.7(6) -lO^^ 


3.3(7) -10^=^ 


(hl6) 


4.1(3) 


0.93(19) 


23.7 


1.10(17) -lO^^ 


8.3(8) -10^=^ 



and for f^,: 



These are deduced from the formulas in ref. |T^, g, with our convention 
jphy steal ^ 131 McV. Bcsides Mr defined in (|l|), the quark mass defined by 
the PCAC relation m^'"^'" also appears in the dimensionless combination 



/i, = m^^^^ro . (5) 

(For the definition of m^*^^*^ see, for instance, the equations (20)-(23) in our 
paper [0].) The dimensionless parameters Btq, frQ and A^tq in (^ and fr^ 
and A^tq in (B) have to be fitted to the data. 



The data from table 3 of ref. and the one-loop ChPT fits are shown 
in figures |^ and |. (The runs with label (c) and (d), which have low statistics 
and hence large statistical errors, are omitted here.) In both figures two fits 
are shown: one taking into account all points and another one where the 
points above Mr = 2 are omitted. As the figures show, the data in the small 
quark mass range clearly show the expected qualitative behaviour with chiral 
logs. The fit parameters have reasonable values, similar to the ones deduced 
in from previous lattice data at larger quark masses |^. For instance, 
the fits with all points correspond to the parameters: Br^ = 8.2, fr^ = 0.27, 
A^tq = 3.5 in formula (||) and frQ = 0.60, A^r^ = 4.3 in formula @. 
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Test of xFT logarithms on 8 xl6 
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Figure 3: Fits of the pseudoscalar meson mass-squared with the one-loop 
ChPT formula. 



The fact that our data in the small quark mass range < rriq < ^rris 
show the expected logarithmic behaviour of chiral perturbation theory is 
quite satisfactory. However, since we are far from the continuum limit and 
we did not yet take into account the Z-factors of multiplicative renormaliza- 
tion, further work is needed to finally deduce the physical values of the fit 
parameters. 



4 Eigenvalue spectra 

The eigenvalue spectrum of the Wilson-Dirac matrix is interesting both phys- 
ically and from the point of view of simulation algorithms. From the algo- 
rithmic point of view the knowledge of low-lying eigenvalues is crucial for the 
optimization of polynomial approximations in TSMB. An interesting ques- 
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Test of xPT logarithms on 8 xl6 
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Figure 4: Fits of the pseudoscalar meson decay constant with the one-loop 
ChPT formula. 



tion is whether there is a statistically significant presence of configurations 
with negative determinant, i. e. with an odd number of negative eigenval- 
ues of the Wilson-Dirac matrix. A significant number of configurations with 
negative determinant would be a serious obstacle for performing numerical 
simulations with and odd number of flavours - as it occurs with m-, d- and 
s-quarks in nature. The results of our collaboration show |T^, IT3] that in 



the investigated range of quark masses no sign problem occurs because the 
statistical weight of configurations with negative determinant is negligible 
even at our smallest quark masses. 
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5 Conclusion and outlook 



Our cost estimates for numerical simulations of QCD with light dynamical 
Wilson-quarks - as a function of the quark mass and volume - allow a rea- 
sonably precise prediction of the computational costs for future simulations. 
Up to now we did not yet perform a systematic investigation of the slow- 
ing down with decreasing lattice spacing, expressed by the power Za in (|l]), 
but previous experience in the supersymmetric Yang-Mills theory [|l^, [T^ at 
small lattice spacings a ~ 0.06 fm and our preliminary results in QCD near 
a ~ 0.17fm are consistent with the expected behaviour Za = 2 [|10|. 

Assuming Za = 2 and starting from the results of runs (h) and {hl6) in 
table 1^ one can deduce the following estimates for obtaining 100 independent 
gauge configurations: 

• On a 24^-48 lattice with lattice spacing a = 0.125 fm at quark mass pa- 
rameter Mr = 0.8 corresponding to m^^ ~ ^m^ and physical extension 
characterized by Lm^ ~ 6 the cost would be C ~ 9 ■ 10^^ — 3 • 10^^ flop, 
depending on whether one considers r^^" or rf^'^'^ as relevant. 

• Similarly, on a 32^ • 64 lattice with lattice spacing a = 0.06 fm at the 
same quark mass and LrriT, ~ 4 the cost would beC~6-10^'^ — 4- 
10^9 flop. 



Compared to most of the estimates in these numbers are lower. Note that 
the costs for parameter tuning and equilibration - which are not negligible 
- are not included in these estimates. Nevertheless, the necessary order of 
magnitude of CPU times for these kinds of calculations is expected to be 
available for the lattice community in forthcoming years. 
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